

AM A A 


AIAA 2003-5495 

LPV Antiwindup Compensation for 
Enhanced Flight Control Performance 

Bei Lu and Fen Wu 

Dept, of Mechanical and Aerospace Engineering 
North Carolina State University 
Raleigh, NC 27695 

SungWan Kim 

Dynamics and Control Branch 
NASA Langley Research Center 
Hampton, VA 23681 -21 99 


AIAA Guidance, Navigation, and Control 

Conference 

August 11-14, 2003/Austin, TX 


For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 
1801 Alexander Bell Drive, Suite 500, Reston, VA 20191-4344 



LPV Antiwindup Compensation for Enhanced 
Flight Control Performance 

Bei Lu and Fen Wu* 

Dept, of Mechanical and Aerospace Engineering 
North Carolina State University 
Raleigh, NC 27695 

SnngWan Kind 
Dynamics and Control Branch 
NASA Langley Research Center 
Hampton, VA 23681-2199 

In this paper, we propose a saturation control scheme for linear parameter- varying (LPV) 
systems from an antiwindup control perspective. The proposed control approach is advantageous 
because it can be thought of as an augmented control algorithm from the existing control system. 

Moreover, the synthesis condition for an antiwindup compensator is formulated as a linear 
matrix inequality (LMI) optimization problem and can be solved efficiently. We have applied the 
LPV antiwindup controller to an F-16 longitudinal autopilot control system design to enhance 
aircraft safety and improve flight quality in a high angle of attack region. 


Introduction 

T HE flight control system of a tactical aircraft has 
different performance goals for low angle of attack 
and high angle of attack regions. For example, in a low 
angle of attack scenario, pilots desire fast and accurate 
responses for maneuver and attitude tracking, whereas 
in a high angle of attack region, the flight control 
emphasis lies in the maintainability of aircraft stabil- 
ity with acceptable flying qualities. The potential of 
high angle of attack flight presents many challenges 
to the control designers. Due to aerodynamic surface 
saturation and control surface effectiveness, unconven- 
tional actuators such as thrust vectoring are suggested 
for aircraft maneuvering at and beyond the stall an- 
gle of attack. However, incorporation of additional 
thrust vectoring hardware could complicate the design 
of flight control laws in the poststall regime. 1 Robust 
multivariable control methods have been recently ap- 
plied to a variety of aircraft models 2 to demonstrate 
their abilities to fly at high angles of attack (see Ref. 
3 and references therein). Besides control law design, 
another major issue of high angle of attack flight is 
control saturation. It is well recognized that actua- 
tor saturation degrades the performance of the flight 
control system and may even lead to instability. The 
destabilizing effects of actuator saturation have been 
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cited as contributing factors in several mishaps in- 
volving high performance aircraft. 4 For this reason, 
various methods of preventing instability due to sat- 
uration have been examined, which include allocating 
control effectors and command scaling and prioritiza- 
tion. 

A popular approach to control saturation is the 
antiwindup method that employs a two-step design 
procedure. The main idea of antiwindup control is 
to design the linear controller ignoring the saturation 
nonlinearities first and then add antiwindup compen- 
sation to minimize adverse effects of the saturation 
on closed-loop performance. Desirable design require- 
ments for antiwindup compensation subject to actua- 
tor saturation are the closed-loop system stability, re- 
covery of the linear design specifications in the absence 
of saturation (linear performance recovery), and the 
smooth degradation of the linear performance in the 
presence of saturation (graceful performance degrada- 
tion). Like other saturation control techniques, the 
antiwindup compensator design often assumes a linear 
time-invariant (LTI) plant and models the saturation 
block as a sector-bounded nonlinearity. Then absolute 
stability conditions (such as Popov, circle theorems) 
are applied for the stability and performance analy- 
sis. 5 

A general framework that unifies a large class of ex- 
isting antiwindup control schemes in terms of two ma- 
trix parameters was proposed in Ref. 6. This frame- 
work is useful for understanding different anti windup 
control schemes and motivates the development of 
systematic procedures for designing antiwindup con- 
trollers that provide guaranteed stability and perfor- 
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mance. Early results in antiwindup control often have 
the drawback of lacking rigorous stability analysis and 
clear exposition of performance objectives. Using an 
extended circle criterion, the synthesis condition of 
static antiwindup controllers is formulated as a linear 
matrix inequality (LMI) problem in Ref. 7. A recent 
study in Ref. 8 has further revealed that antiwindup 
control for stable open-loop LTI systems can be solved 
globally as an LMI problem with the order of anti- 
windup compensator no more than the plant’s order. 
Alternatively, the Popov stability condition has also 
been applied to the antiwindup compensator design 
problem. 9 However, the synthesis condition of the sat- 
uration controller is given in coupled Riccati equations, 
which is difficult to solve for its optimal solution. Most 
previous antiwindup compensator designs are only ap- 
plicable to open-loop stable LTI systems, limiting their 
usefulness for practical problems. When the system is 
nonlinear and open-loop unstable, the control synthe- 
sis problem becomes very difficult to solve, therefore, 
global stabilization cannot be achieved. 10, 11 However, 
in many control systems including flight control sys- 
tems, the system dynamics are inherently nonlinear 
and their linearizations are strictly unstable. 

The motivation for this research is twofold. First, 
the antiwindup control scheme for LTI plants in Ref. 
12 is generalized to linear parameter-varying (LPV) 
systems. This generalization is very important be- 
cause of the relevance of LPV systems to nonlinear 
systems. In fact, the LPV model can be thought of 
as a group of local descriptions of nonlinear dynam- 
ics. The antiwindup compensation augments existing 
control systems by maintaining stability and recov- 
ering control performance when actuator saturated. 
Second, saturation control for aircraft under large ma- 
neuver operations is critical due to safety concerns. 
The proposed antiwindup compensation can be devel- 
oped by augmenting existing flight control algorithms. 
The consequence is enhanced reliability and an ex- 
panded flight envelope. In particular, we demonstrate 
adequate flight control performance using LPV anti- 
windup control in a high angle of attack scenario. 

The notation in this paper is standard. R stands for 
the set of real numbers and R + for the non-negative 
real numbers. R mx ™ is the set of real m x n ma- 
trices. The transpose of a real matrix M is denoted 
by M T . Ker(M) is used to denote the orthogonal 
complement of M. A block diagonal matrix with 
submatrices Xi,X 2 , . . . ,X p in its diagonal is denoted 
by diag {Xi, X 2 , . . . ,X P }. We use S” xn to denote 
the real symmetric n x n matrices and S” xn to de- 
note positive definite matrices. If M £ S" xn , then 
M > 0 (M > 0) indicates that M is positive definite 
(positive semidefinite) and M < 0 [M < 0) denotes 
a negative definite (negative semidefinite) matrix. If 
a,b £ R, then sect [a, b] denotes the conic sector de- 
fined by {(q,p) : (p — aq){p — bq) < 0}. For x £ R”, 


its norm is defined as ||x|| := (x T x)i. The space of 
square integrable functions is denoted by C 2 , that is, 

for any u G C 2 , ||tt||2 := [/ 0 °° u T (t)u(t,)dt ] 2 is finite. 

LPV Antiwindup Control Synthesis 

Antiwindup compensation is to modify nominal con- 
trollers so that if the signal from the controller is 
different from what enters the plant, corrective feed- 
back action is employed to reduce the discrepancy. 
Because it is impossible to provide a global stabiliz- 
ing solution to the antiwindup control problem when 
the open-loop plant is unstable, one often needs to de- 
termine regional stability for saturation control and 
to design controller gain in the guaranteed stability 
region. 10,11 In Ref. 8, a sector-bounded input non- 
linearity, sect [0,1], was considered for the open-loop 
stable plant and is not applicable to exponentially un- 
stable systems. However, the derived performance and 
stability properties can be improved when the input 
nonlinearity is restricted to a smaller sector region. As 
a result, this modification leads to regional stability of 
the antiwindup compensated system and extends the 
antiwindup control technique to exponentially unsta- 
ble open-loop systems. 

Consider an LPV plant P p described by 


Xp 


Ap(p) Bpi(p) B p 2 {p) 


x p 

e 

= 

Cpilp) D p hIp) D pl 2 {p) 


d 

y _ 


_Cpi{p) Dp 2 ilp) D p 22 (p)_ 


_cr(u)_ 


where the plant state x p £ R np . y £ R' l » is the 
measurement for control, and cr{u ) £ R’ 1 " is the satu- 
rated control input, e G R ne is the controlled output 
and d £ R n<1 is the disturbance input. It is assumed 
that the vector-valued parameter p evolves continu- 
ously over time and its range is limited to a compact 
subset V C R s . In addition, its time derivative 
is assume to be bounded and satisfy the constraint 
Hi < Pi < V , i = 1, 2, . . . , s. For notational purposes, 
denote V = {v : za < v < P*, i = 1,2,..., s}, where V 
is a given convex polytope in R s that contains the ori- 
gin. Given the sets V and V, the parameter u - variation 
set is defined as 

^ = {pG C 1 (R+,R S ) : p(t) £ V , p{t) £ V, Vt > 0} 

So specifies the set of all allowable parameter tra- 
jectories. In the interests of notational compactness, 
the parameter dependence will not always be shown 
from now on. 

All matrix valued state-space data are continuous 
and have appropriate dimensions. For simplicity, we 
assume that 

(Al) (A p , B p2 ,C p2 ) triple is parameter-dependent 
stabilizable and detectable for all p. The LPV 
system (1) is said to be parameter-dependent sta- 
bilizable if there exists continuously differentiable 
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matrix functions X(p) = X T (p) > 0 and F(p), 
such that for (p, p) € V x V 


(A p + B P 2F) T X + X(A p + B P 2F) + 

»= 1 


. dX 
Pijr~ < 
O Pi 


o 


(A2) The matrices \B p2 an d \C P 2 -D p2 i] 

have full row rank. 

(A3) D p 22 = 0. 


are used to condition the nominal controller when the 
control input is saturated. 

Our objective is to design an LPV antiwindup com- 
pensator AW P such that the adverse effect of input sat- 
urations are minimized in terms of induced £2 norm. 
The antiwindup compensator is in the form of 


Aaw iPi p) 
Caw iPi p) 


B a w {p) 

D a w ip) 



%aw 


q 


( 3 ) 


The actuator nonlinearity under consideration is a 
piecewise-linear saturation 


f Ui |lii|<<"“ 

\ sgn(wj)u™ ax \ui\ > u“ ax 


for i = 1, 2, . . . , n u . The antiwindup control structure 
is shown in Fig. 1(a). 



(a) LPV antiwindup controller structure 



(b) Equivalent transformation 

Fig. 1 Nonlinear saturation control diagram. 

Following the standard antiwindup procedure, a 
nominal LPV controller K p is designed first by ig- 
noring the input nonlinearity. Different linear control 
design techniques can be employed to achieve this goal. 
A systematic way to do this is through LPV control 
theory. 13-15 Due to the assumption (Al), the nomi- 
nal controller K p is capable of stabilizing the open-loop 
system when no input saturation exists, and its design 
will determine the nominal performance of the closed- 
loop LPV system. We assume that such a controller 
is given by 


Xk 


Akip,p) 

Bkip) 

Xk 

1 

Vl 

u 


. Ckip) 

Dkip) 

_y _ 


y>2_ 


where Xk G R rafe . The variables v \ , V 2 are the auxiliary 
inputs provided by an antiwindup compensator. They 


where the state x aw € R n ™; the size of the compen- 
sator state is determined later on. Such antiwindup 
compensation schemes provide a computationally effi- 
cient technique for the “retro-fitting” existing uncon- 
strained controller K p to account for input nonlinear- 
ities, thereby eliminating controller windup problems 
for input saturated nonlinear systems. 

The LPV antiwindup control diagram in Fig. 1(a) 
can be transformed to its equivalent form by substitut- 
ing the actuator saturation with a deadband nonlinear- 
ity Aj = 1 - and A = diag{Ai, A 2 , . . . , A„ u }, 
as shown in Fig. 1(b). This allows recasting of 
the compensator design problem into a robust control 
paradigm. 



Fig. 2 Restricted range input saturation. 


It can be seen that the deadband uncertainty A, 
resides in the conic sector [0,1]. The magnitude of 
the control input signal u t is restricted to be less 
than ^ j M » x by reducing the nonlinearity A, to 

sect [0, A;] with 0 < k < 1. Consequently, we have 
A G sect[0, kl\. The resulting stability notion then 
becomes regional stability. However, this restriction 
of uncertainty A extends the antiwindup scheme to 
open-loop exponentially unstable systems. 

Let the system G p be the interconnection of the 
open- loop system P p and the nominal controller K p , 
but exclude the antiwindup compensator. Then its 
dynamic equation is 


X 


u 


e 


Q_ 



A(p, p) 
Co(p) 
Ci(p) 
0 


q = A u 


Bo{p) Bi{p) 

Dooip) An(p) 
Dioip) Diiip) 
I 0 


Bi(p) 
D 02 ip) 
Dl2ip) 
0 


1 x 

q 

d 

Vl 

J v 2 


( 4 ) 

( 5 ) 
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where x £ R" with n = n p + n k and A £ sect[ 0 , kl], 
A = 


A p T B p 2 D k C p2 B p 2 C k 
BkC p2 A k 


Bn = 


B 2 = 


-B 

o' 


P2 


B i = 


B p i + B p2 D k D p2 1 
B k D p 21 


0 B p 2 

/ 0 


Co = [D k C p 2 Ck] 


Ci = \C P i + D p \ 2 DkC p2 D p i 2 Ck\ 

Boo = 0 , D 0 i = D k D p 2i 

-£>02 = [0 /] , -Dio = ~D p 12 

-Du = Dpn + D p i 2 DkD p2 i, D\ 2 = [0 -D p i2] 

Note that the state-space data have linear dependency 
of parameter rate p. 

Denote = [ir T x^ u ,] . Then the final closed-loop 
system T p = TAG P , AW p ) is described by 


As expected, the closed-loop state-space data depend 
on the antiwindup compensator gain in affine form. 
The following theorem provides a synthesis condition 
for the antiwindup compensator. 

Theorem 1 (Synthesis Condition for LPV Anti- 
windup Compensator) 

Given a scalar 0 < k < 1 , the LPV open- 

loop system P p with a parameter- dependent stabiliz- 
ing nominal controller K p , if there exists a pair of 
positive- definite matrix functions Rn(p) £ s" pX " p , 
S(p) £ S” x ™ and a diagonal matrix function V(p) = 
diag {ui(p), • • • , v„ u (p)} > 0 satisfying 


RuAl + A p Rii - ^ 1 b p 2 VB ^ 2 - 


Xd 


Acl(P: P) 

Bo, dip) 

Bi,dip) 


x c i 

U 

= 

Co, clip iP) 

Doo.ciip) 

Doi,dip) 


q 

e 


Cl, clip, p) 

Dio, dip) 

Du, dip) 


d 






( 


q = A u. 


( 7 ) 


Gathering all parameters of the antiwindup compen- 
sator into the single variable 


, iipi.. 1± fe2 -Op2 v J -> p2 

CpiRii - ^±D p i 2 VB, 
B T p i 

RuC^i - 2 -^B p 2 VDT 12 

Jpl2v -^pl2 

■'pll 

S r ) <? 

SA + A T S + Yp t — SB 


2=1 


dpi 


T 

p 2 


B 


- 'yin. - 2 S Y^D v i 2 VDl 

dT 


pi 


D 


p li 

- 7 -fn d 


<0 


( 9 ) 


BfS 

Ci 


2=1 


dpi 


“7 J n d 
Du 


Cf 

D T n 

“7 In t 


< 0 



r 

Rn 

[ln P A]' 

0(p,p) = 

A R 

-ci aw J -'aw 

Caw Daw 


In p 

0 

s 


> 0 


(10) 

( 11 ) 


and introducing the shorthands: 

•A(p,p) = 


A 0 

0 0 


for all (p,p) £ V x V, then there exists an n p -th 
order LPV antiwindup compensator AW P to stabilize 
the closed-loop system quadratically and have the per- 
formance || e|| 2 < 7IMI2 for all p(-) £ Tf, when the 


Mp) = 


B 0 

0 


Bi(p) = 


Bi 

0 


condition \uf\ < (izi) u i 


i = 1, 2, . . . , n u holds. 


C 0 {p) = [Co 0] , Ci(p) = [Ci 0] 


vlip) 

Vlip) 

n f(p) 


0 b 2 

1 0 


0 -D, 


02 


[0 

I 

0 

0 ' 

0 
1 

0 


0 


0 -D12 

[Qi 0.2 Q3] 


1 1. closed-loop state-space data related to the inter- 
connected system G p and the antiwindup compensator 
AW P is as follows 


Proof: Consider a Lyapunov function V (x) = 

x T X c i{p)x for the closed-loop system T p , then a suffi- 
cient condition for the stability and performance can 
be established from the inequality 

V + —e T e — r yd T d + 2 q T W ( ku — q) < 0 

7 

using 5 -procedure. Note that it is equivalent to 

dX cl 


A c i 

Bo, ci 

Bi, c i 


Co, cl 

Doo,ci 

Doi,ci 


Cl, d 

Dio, d 

Dn,d_ 


'A 

Bo Bi~ 


A 

> T 1 

1 

Co 

D 00 D01 

+ 

1 

>T 

2 

Ci 

Dio Dn_ 


l- 

>T 

3 


AhX c i 


XciAd 




B( cl X cl + kWCo.d 


e [Qi Q2 q 3 ] ■ ( 8 ) 


BldXci 
Cl, cl 




* 

* 

* 


kiWDoo.d + Dl o cl W) - 2 W 

* 

* 

< 0 

kD^ cl W 

- 7 1 

* 

Dio, d 

Du, d 

-7/ _ 

(12) 
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where * in the upper triangle means symmetric. Us- 
ing the expressions in (8), the inequality (12) can be 
rewritten as follows 

*(p,p)+v T (p)e(p,p)Q+Q T e T (p,p)v(p)< o (13) 

with 


Mq^JVq < 0 are the same as conditions (9) and (10), 
respectively. 

Given the definition for matrices X c i and X~ t , the 
coupling condition between R and S would be 

R I 

. „ > 0 and rank(A — A -1 ) < n aw 


A T X cl + X cl A+Y j p i ^ 

* = B%X cl + kWC 0 

BfX cl 

Ci 

* * * 

k(WD 00 + D£ 0 W) -2 W * * 

kDj y[ W -7 1 * 

Dio -Dii — 7 1 

V = \PxX cl kV 2 W 0 V 3 ] 

Q = [Qi Q 2 Q 3 0] 

Partition the matrix X c i compatibly to the states of 
an interconnected system G p and antiwindup compen- 
sator AW P as n = n p + nk and n aw , and let 


for all p £ V. Because only the (1,1) element of R 
matrix is constrained in the LMIs (9)— (11), it is al- 
ways possible to augment matrix An to R in satisfying 
the above coupling condition. For example, one may 
choose 


A 11 

[/ 0] A - 1 

o' 

I 


[0 i] s- 1 

I 

0 


[0 I] S - 5 1 

o' 

I 



The resulting R matrix is positive definite because of 
the condition (11). Also, R — A -1 > 0 is satisfied for 
the selected A matrix. The rank condition is trivially 
satisfied if one chooses n aw = n p . So we obtain 
the desired synthesis condition for the antiwindup 
compensator. Q.E.D. 


Xd(p) 

Ki\p) 


S(p) 

N(p) 



_N T { P ) 

7 



■ R{p) 

M{p) 



M T (p ) 

7 



An (p) 

Al2(p) 

M(p ) 

R-Up) 

A 22 (p) 

M 

Hp) 


7 


where MN T = I — RS , and “?” means the elements 
“we don’t care”. According to elimination lemma, 16 
the inequality (13) is equivalent to 

Af£(p)^{p,p)Af v (p) < 0 Uq^(p,p)Uq < 0 (14) 


where A fp and Mq are the null spaces of matrices V 
and Q, which are 


M-p = diag 



w-\ij 


' I 

0 

0 

0 ' 

0 

I 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I 

0 

0 

0 

0 

/ 


I 

0 

0 ' 

0 

0 

0 

0 

0 

0 


— D 1 

^p 12 

0 

0 

0 

I 

0 

/ 

0 _ 


Denote V(p) = W 1 (p). Through lengthy algebraic 
manipulations, it can be shown that Afp ftfp < 0 and 


Remark 1 The minimum order of an antiwindup 
compensator that can always he synthesized by a con- 
vex condition is n p . 8 However, if we enforce R = S ~ 1 , 
then it is possible to solve a modified LMI optimization 
problem from (9)-(ll) 


DnAj + A p R n - 2 -^B p2 VB T v2 




dRu 

dpi 


C pl R n - 2 -^D pl2 VBf 2 

BL 


R 11 C £ - 2 ~^B p2 VDl 12 B p i 

-7 / - ^^#D pl2 UDj 12 D p ii 

Dl ii -7 1 


ar+ra t -y>P 

'H dpi 

B\ 

RCf 


L _L 

BT 

- 7 / 

Bu 


Cl A 

D n 

-7 1 

II 

'533 

R 12 } > 0 
^ 22 




< 0 (15) 


<0 (16) 
(17) 


This results in a static LPV antiwindup compensator. 
However, the modified LMIs do not necessarily have a 
solution. 


Recall that the LPV antiwindup compensator syn- 
thesis problem is originally formulated as a robust 
control problem. This usually leads to a noncon- 
vex solvability condition as bilinear matrix inequalities 
(BMIs). However, due to the special structure of the 
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antiwindup controller, the resulting synthesis condi- 
tion is convex in terms of matrix variables Rn,S, 
and V. In fact, the solvability condition for the LPV 
antiwindup compensator is provided as an infinite- 
dimensional LMI optimization problem, for which an 
efficient numerical algorithm exists to solve it approx- 
imately. 17 This can be done by parameterizing the 
matrix variables using a finite set of scalar basis func- 
tions as 


Nf 

Rn(p) = 5>(p)Jlii t « 

i=l 

N h 

V (P) = ^2,h k {p)V k 

k = 1 


Ng 

S(P) = ^9j{p)Sj 
i=i 


where {fi(p)}?J x , {gj(P)}^=i and { h i(p)} k =i are user - 
specified scalar basis functions. Ru,i,Sj, and V k are 
new optimization variables to be determined. After 
such a parameterization, the LPV synthesis conditions 
can be solved using a gridding method over the param- 
eter space. 

After solving Rn,S, and V matrix functions, the 
LPV antiwindup compensator gain can be determined 
by either solving the feasibility problem from the 
closed-loop LMI (13) or the antiwindup compensator 
can be constructed explicitly as shown in the follow- 
ing theorem. The explicit construction approach is 
advantageous because it avoids possible numerical ill- 
condition when solving the above feasibility LMI prob- 
lem. Moreover, it connects the antiwindup controller 
directly to the plant and nominal controller gains. 


Theorem 2 (LPV Antiwindup Compensator Con- 
struction) 

Given the solutions Rn,S, 7 and V = W~ x of the 
LMIs (9)-(ll), let MN t = I n - RS with M,N £ 
graxnp an( } £T _ | ~j n ^ 0 ] , then one n p -th order 
LP V antiwindup compensator AW P can be constructed 
through the following scheme: 


1. Compute a feasible D aw (p) G R n “ xn “ such that 


n u x rip 


0 

In u 0 O' 



0 

-n 

0 



Bl 


0„ 


S + kWC 0 
BfS 

Ci 


0 

kW T 0 Dj 12 ■ 

kW 


0 

n 

Dp\2 



C a 


B T p2 


[Bl + kWC 0 RE + Df w Bf 2 
BlE 
Ci RE 


and the matrix A aw (p, p) £ R nx ™p as 
A aw = -p^-RE - p^M T E - A t E 

dp dp 

- { SB 0 + B aw + kC^W SB 1 Cf] n- 1 

'’(Bg’ + kWC 0 R)E + b T aw B T p2 + kWC av 

BfE 

C\RE + D p i 2 C aw 


3. Convert the transformed antiwindup compensator 
gain to its original state-space data by 


'A U 

aw ±J aw 


'N SB 2 

Caw Daw 


,° [0 In u \_ 


A* 


Ba 


C F) 
p-'aw Jy aw . 

m t e 0 

0 In.. 


SARE 0 

0 0„..v. 


Proof: The derivation of the antiwindup controller 
formula basically follows the procedure outlined in Ref. 
12. Define 


-k(WD aw + D T aw W) + 2W * * 

'kB'nW ^ rind * 

— (Dio + Dpi 2 L) aw ) —D 11 7 I„. e 


> 0 


for all p £ V. Denote the above inequality as 

n (p). 


z 1 


I RE 

0 M t E ’ 


Z2 


S E 

n t 0 


Then it can be shown that X c iZi = Z 2 . Also we have 
the following congruent transformation: 

ZlX cl Zi = 


S E 
E t E t RE 


2. Compute the least square solutions of the following 
linear equations for B aw (p) £ R nxn “, C aw (p) £ 



dS_ 

dt 


E T (R d l + M^r) 


(R d l + M^r) T E 

-E T ^E 
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ZfX c iA c iZi ZfX cl B 0 , c i ZfXciBi’Ci 

Co,ciZi Dqo, c i -Doi,cZ 

Ci ,clZ i -Dio.cZ Dii, ci 


' SA 

0 

SB 0 

SBi 

E t A 

e t are 

E T B 0 

E t Bi 

C 0 

C 0 RE 

Dqo 

^01 

Ci 

Ci RE 

Dio 

Du 


In 0 


' A 

-P^-aw 

Baw 


' 0 I np 

0 

0 ' 

Caw 

Baw_ 


0 0 

I 

0 


0 -Dp 12 

where 


The (1, 1) element of the above matrix inequality 
corresponds to LMI (10) after elimination of the 
variables B aw and D aw . It can also be shown that the 
(2, 2) element is equivalent to LMI (9) by eliminating 
CawiD aw . Q.E.D. 

The explicit antiwindup construction scheme can 
also be applied to open-loop stable systems. How- 
ever, because the synthesis condition for an open-loop 
stable plant does not involve matrix IT, 8 we need to 
solve the feasible D aw and W matrices altogether at 
the first step. The remaining steps are the same by 
setting k = 1. 


~ A 

-flaw 

Baw 


\SARE 0] 

r 1 

aw 

B aw _ 


o 1 

o 
1 


'N SB 2 ' 

'A 

aw 

Baw 

\m t e ol 

_0 [0 /]_ 

Caw 

Baw 

o 
1 


Multiply diag {Zj , /, /, /} from the left side and its 
conjugate transpose from the right side of eq. (13), 
and we get 


' Hu 

* 

* 

H 2 1 

h 22 

* 

L\ 

L 2 

-n 


< 0 


(18) 


with the shorthand notation 


dS 


H 1 i = SA + A T S + Y j p. 

U d P 


dS 


dN 1 


H 2 i — E A + A aw + E ( R — + M ) 


dt 


dt 


dR, 


H 22 = E t {AR + RA t - Vp t -)£ 

V Pi 


+ B p2 C a w + Ca W Bp 2 


L i — 


BZS + B t 


aw 1 

B^S 

Ci 


kWC 0 


L‘2 — 


\{Bl + kWC 0 R)E + D T aw B T V 2 + kWC a 
BfE 

GiRE + D p i 2 C aw 
By Schur complement, it is equivalent to 

< 0 (19) 


Hu + Li n ~ L Li 

21 T JJ 2 - 


H 2 i + 1 Li h 22 


lIa~ 1 l 2 


Clearly, the lower (3 x 3) matrix of the inequality 
(18) must be negative definite, this determines the fea- 
sible D aw . Let the (2, 1) element be equal to zeros and 
one can solve for A aw . This also leads to decoupled 
LMIs from the inequality (19). Then B aw ,C aw terms 
can be solved from the (1, 1) and (2, 2) elements of the 
decoupled inequality (19). Note that both inequalities 
have regular solutions. 18 


Saturation Control for Flight Dynamics 

In this section, the proposed LPV antiwindup con- 
trol synthesis technique is applied to flight dynamics. 
The system to be controlled is the longitudinal F-16 
aircraft model based on NASA Langley Research Cen- 
ter (LaRC) wind tunnel tests, 19 which is described by 
Stevens and Lewis in great detail. 20 The full nonlin- 
ear longitudinal model of an F-16 aircraft is given as 
follows: 


• qSc . , cos a 

V = - — [C xq cos a + C zq sin a\ q + 1 


2 mV 


qS 


— g sin(6 l — a) + — \C Z C os a + C z sin a] (20) 

TO 


qSc 

1 + 2mC 2 ^ Czq C ° S a ~ Cxq Sin a ) 


q-T 


,sin a 

mV 


+ ^ cos(0 — a) + [ C z cos a — C x sin a ] (21) 


qSc , qSc 

q — 2 ^ y [cCmq- + XC zq \ q + —j— 


= q 


C m + ~C Z 

C 


(22) 

(23) 


The aerodynamic coefficients are provided as lookup 
tables from NASA LaRC wind tunnel tests on a scale 
model of F-16 aircraft. The data apply to the speed 
range up to about Mach 0.6, and cover a very wide 
range of angle of attack (—20° < a < 90°). How- 
ever, conventional aerodynamic math models for use 
in aircraft simulation or flight control design have be- 
come increasingly deficient in the poststall region. So, 
the investigation on the robust control of aircraft in 
the high angle of attack region with uncertain aero- 
dynamic coefficients is one of the challenging research 
topics. 

The F-16 is powered by an afterburning turbofan 
jet engine, which produces a thrust force in the x axis. 
The NASA data include a model of the engine in which 
the thrust response is modeled with a first-order lag; 
the lag time constant is a function of the actual engine 
power level and the commanded power. The throttle 
position is related to the commanded power level. For 
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convenience, the actual power level is also considered 
as a state variable in longitudinal dynamics. 

The nonlinear model is then linearized at trim con- 
ditions to generate state-space models for design. The 
state and input variables are as follows: 


x = < 


V 

a 

Q 

e 

pow 


true airspeed (ft/s) 
angle of attack (rad) 
pitch rate (rad/s) 
pitch angle (rad) 
actual power level (0-100) 


u = 



throttle position (0-1) 
elevator angle (deg) (— 25°-25°) 


By slight abuse of notations, the above listed variables 
are perturbations from their equilibrium states. In ad- 
dition, V, g, and flight path angle 7 are selected as 
outputs. The operating envelope of interest covers air- 
craft speeds between 160 ft/s and 200 ft/s and angles of 
attack between 20° and 45°. These two variables are 
used as scheduling parameters in the LPV modeling 
of F-16 longitudinal dynamics. Fig. 3 shows the flight 
trim points for the control design. The corresponding 
dynamic pressure of the designed range is between 20 
psf and 50 psf. This portion of flight envelope is cho- 
sen because the moderately high angle of attack and 
the low dynamic pressure cause aerodynamic control 

surface saturation. 

18000 

15000 
12000 

& 

■§ 9000 
6000 
3000 



Fig. 3 Flight trim points. 

The design objective of the nominal LPV controller 
in this research is to track the flight path angle com- 
mand with the tracking error about 1.25% of the com- 
mand in the steady state. This kind of problem is 
formulated as a model-following problem, 21,22 where 
the ideal model to be followed is chosen to be a second- 
order filter based on desired flying qualities 

2 

7ideal ^ideal (24) 

7cmd d 2 + 2/ it l ea lU'i ( | fta | S + uf deal 

The implicit model-following framework allows for di- 
rect incorporation of flying quality specifications into 
the control design. A block diagram of the system in- 
terconnection structure for synthesizing the nominal 
controller is shown in Fig. 4. 



Fig. 4 Open-loop interconnection for LPV control 
design. 

The weighting functions are chosen as 


W p = 


80(s/5 + l) 
s/0.05 + 1 


W n = diag{0.8, 0.6, 0.1} 


lU = diag{l,10,i,+} 

T „ _ 2.25 

ideal s 2 + 2.4s + 2.25 


The throttle and elevator actuator dynamics are mod- 
eled as first order filters. 

8th = _±_ Je _ = 20 

<5™ d s + 5 (5/ md s + 20 

Both the position and the rate of control inputs are 
fed into W u to penalize the control effort. Therefore, 
the system matrix of W ac t is derived as follows. 


W'act = diag 



The nominal controller is designed by formulating 
an LPV synthesis problem, which can be solved us- 
ing either a single or parameter-dependent quadratic 
Lyapunov function over all gridding points in the two 
dimensional parameter space. 13,14 In this work, the 
single quadratic Lyapunov function (SQLF) approach 
is chosen to reduce the computation time in the nom- 
inal LPV control synthesis phase. The performance 
obtained through the SQLF approach is 7 n0 m = 8.0. 

The nonlinear F-16 longitudinal dynamics with the 
nominal LPV controller K p is simulated, and the mag- 
nitude limits of actuators are enforced during the sim- 
ulation. One designed flight condition is used as the 
test point, which is a trimmed flight at V = 160 ft/s 
and a — 35°. Figures 5 and 6 show the nonlinear 
time responses of the aircraft model to the flight path 
doublet input with different magnitudes. 

It is observed that when flight path doublet input 
is equal to ±1°, the positions of throttle and eleva- 
tor are not saturated. The simulation result of the 
flight path shows that the measurement matches the 
ideal response very well. For comparison, the nonlin- 
ear time response for a more drastic ±2° flight path 
doublet input is shown in Fig. 6. It should require 
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more control actions of the actuators to achieve the 
performance objective. Unfortunately, the elevator is 
saturated severely and thus makes the system unsta- 
ble. 

The proposed LPV antiwindup scheme is applied to 
the control of the F-16 aircraft, and we also choose 
the SQLF approach to perform the LPV synthesis. 
A series of LPV antiwindup compensators AW P are 
designed by gradually decreasing k value from 1 to 
0.90. Table 1 shows the corresponding performance 
level 7. Because one of the gridding points is unstable, 
the synthesis condition for LPV antiwindup control 
is infeasible for k = 1. Compared with the nom- 
inal performance, the worst performance level when 
k =1 — 10 -5 indicates the strong adverse effect of sat- 
uration nonlinearity. When k = 0.90, the antiwindup 
compensator almost recovers the nominal closed-loop 
performance 7 nom . 

Table 1 11 -^ performance level vs. sector range 

[0 ,kl}. 


Sector range [0, kl] 

Uoo performance 7 

1 - 1CT 5 

1596.38 

1 - 10 -4 

258.30 

1 - 1CT 3 

32.16 

0.99 

8.41 

0.90 

8.00 


Now, a ±2° flight path command input is applied 
to the F-16 aircraft system with the antiwindup com- 
pensator. The time response of the nonlinear system 
at the same trimmed condition is tested. As shown in 
Fig. 7, the antiwindup compensator not only stabilizes 
the system when saturation occurs, but also provides 
fairly good recovery of tracking performance. 

Conclusion 

We have developed an antiwindup control scheme 
for LPV systems and applied the proposed saturation 
control scheme to F-16 aircraft. Due to the special 
structure of the antiwindup control scheme, the LPV 
antiwindup control synthesis condition is solvable by 
LMI optimizations. The extension of the antiwindup 
control idea to LPV systems provides a practical ap- 
proach for nonlinear flight dynamics in the presence 
of actuator saturation. Saturation control is par- 
ticularly important to near stall and poststall flight 
conditions. By augmenting the nominal longitudinal 
autopilot with antiwindup compensator, it was shown 
through nonlinear simulation that the F-16 aircraft 
maintains its stability and has adequate control perfor- 
mance in a high angle of attack. Applying antiwindup 
control scheme to flight control is promising because no 
additional actuator is needed to compensator control 
authority. In fact, the implementation of antiwindup 
controllers could be done by modifying flight control 
software. 
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(a) Flight path angle 



(b) Angle of attack 



(c) Elevator angle 



(d) Throttle position 


(a) Flight path angle 



(b) Angle of attack 



(c) Elevator angle 



(d) Throttle position 


Fig. 5 Nonlinear simulation of F-16 dynamics to 
±1° flight path doublet input. 


Fig. 6 Nonlinear simulation of F-16 dynamics to 
±2° flight path doublet input. 
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(a) Flight path angle 



(b) Angle of attack 



(c) Elevator angle 



(d) Throttle position 


Fig. 7 Saturation control of F-16 dynamics using 
LPV antiwindup compensator. 
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